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Abstract 

We consider universal statistical properties of systems that are characterized by phase states 
with macroscopic degeneracy of the ground state. A possible topological order in such systems is 
described by non-linear discrete equations. We focus on the discrete equations which take place in 
the case of generalized exclusion principle statistics. We show that their exact solutions are quantum 
dimensions of the irreducible representations of certain quantum group. These solutions provide an 
example of the point where the generalized exclusion principle statistics and braid statistics meet 
each otiicr. Wc propose a procedure to construct the quantum dimer models by means of projection 
of the knotted field configurations that involved braiding features of one-dimensional topology. 



PACS numbers: 71.10.-w, 71.10.Pm, 71.10.Fd 



I. INTRODUCTION 



The universal behaviour of low-dimensional strongly correlated systems at low temper- 
atures is determined to a great extent by the topology of the manifolds, where the ground 
state and low-lying excitations are determined. In strongly correlated electron liquids with 
high degree of degeneracy of the ground state such a manifold is presented by a collection 
of strings in the form of an arbitrary tangle of knotted and linked filaments. The processes 
of fusion and decay of the strings give rise to modifications of the tangle and, consequently, 
to the universal character of the quantum criticality in phase states with topological order. 
Being a result of detailed study of electron hquids in the states with fractional Hall effect, 
this conclusion was supported recently by the results, obtained during studying dynamics 
of spinP^El g^j^j charg^^ degrees of freedom in other low-dimensional electron systems. 

It is well known that statistics of excitations in (1 + 1)- and (2 + l)-dimensional systems is 
connected with the braid group. Quantum states in such systems are classified by irreducible 
representations of the braid group instead of the even (for bosons) or odd (for fermions) 
irreducible representations of the permutation group as (3+l)-dimensional systems. One- 
dimensional irreducible representations of the braid group correspond to Abehan anyon 
states, while multi-dimensional irreducible representations describe the non-Abelian states of 
anyons. Statistics of anyon excitations is called either braid statistics or fractional statistics, 
because it leads to the conclusion that there exist particles with a fractional charge and spin. 
In the long-wavelength limit the description of non-Abelian anyons is based on the effective 
action of the topological field theor}^, which contains the Chern-Simons term. An important 
particular case is the quantum group SU{N)k case, where the integer k is the coefficient of 
the Chern-Simons action, or a level in the Wess-Zumino-Witten-Novikov theory. 

We can also employ another approach'^, based on the generalized exclusion principle. 
Studies of the distribution function using this approaclP^ninilini j^^ve shown that this 
method is equivalent to the thermodynamic Bethe ansat^. The equation, which determines 
the minimum of free energy, has the form of the Hirota equation!^. It is well known from the 
theory of nonlinear equations, that this discrete equation in the continuous hmilP^ yields the 
known integrable hierarchies of nonlinear equations. An important feature of the derivation 
of discrete equations in the theory with the generalized exclusion principle is the absence of 
any reference to the dimensionality of the space, unhke braid statistics. Formally generahzed 
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exclusion principle statistic^ may take place not only in low-dimensional systems. The 
reason of the emergence of the solutions, which exist only in low-dimensional systems, is 
that in the limit of large values of momentum, discrete equations in the universal sectoif^ 
of exclusion statistics encode actually the particle fusion rules, which take place in the 
conformal field theory. Going over to the limit /c ^ 1 in discrete equations of motion, the 
specifity of low-dimensional situation vanishes. 

The features of considered statistics are summed up by the theory of tensor categorie^^S^ll. 
It unifies consistently the processes of braiding and fusion of string manifolds, which are 
images of quasiparticle world fines. In the case of low-dimensional systems this occurs in the 
theory of modular tensor categorie^ j^^'^^ l The theory of symmetric tensor categorie: j^^'^" ' ^^' is 
applicable to the (3 + 1)D systems. In the last case, restrictions on solutions of equations 
of motion are so strong, that aU anyon states are excluded, and only bosons (a = 2tt) and 
fermions (a = vr) with the interchange phase a remain out of anyon states. 

To solve the problems of the theory of strongly correlated systems on a lattice, it is 
often convenient to use the theory of the braid group representation or the Temperley- 
Lieb algebra (TLA) representation^ with a special value of the TLA parameter. In the 
continuous limit, we can also employ the effective Chern-Simons actionP^'^I! in particular, 
to classify the hierarchy of phase states with the aid of (3 + 0)-dimensional spinor Ginzburg- 
Landau functionaP^, it is convenient to use the Hopf invariantPSEMHl^ which is the {3 + 0)D 
analog of the Chern-Simons term. 

The construction of a lattice model out from the continuous theory (even inheriting its 
essential properties) is evidently an ambiguous procedure. Some intuitive insight of how this 
can be done, based on the mentioned properties, may include the foUowing consideration. 
It is well known that the Chern-Simons term in the action of (2 + l)-dimensional systems 
encodes the invariant description of fiuctuating Aharonov-Bohm vortices. The action of 
the doubled Chern-Simons theory (for example SU{N)k x SU{N)p^^^^^ which refers to 
the systems, where T and P invariances are not brocken, includes pairs of Aharonov-Bohm 
vortices with the opposite chiralities. It is natural to suppose that neutral pairs of Aharonov- 
Bohm vortices in such systems are plane slices of a string loop (from the three-dimensional 
point of view). A pair of Aharonov-Bohm vortices are portions of the loop cutted by a plane. 
This means on the whole that a small loop with the scale of the order of the lattice constant 
induces a dimerized configuration of currents when it is projected on the plane. Such a 
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projected loop, or equivalently the dimer configuration, can be a building block for the 
formation of self-organized mesoscale structures in the form of nets. The increased interest 
to quantum dimer distribution^^^'^SEQ] jg connected not only with the theory of resonance 
valence bonds or with the support originating from the experience of the exact solvable 
model^, but it is also motivated by recent result^^Ml the field of non-Abelian gauge 
theory. 

The goal of this paper is two-fold. In the second section we will discuss the opportuni- 
ties, which appear due to mapping essential configurations of the (3+0)-dimensional spinor 
Ginzburg-Landau model into lattice (2 + 1)D dimer configurations. We will show that two- 
dimensional dimer field configurations, distributed on two sub-lattices, may be obtained by 
projecting three-dimensional current configurations of the so-called toroid phase state. In 
the third section we will consider the solutions of nonhnear discrete equations of the ther- 
modynamic Bethe ansatz and will show their relation with the characteristics arising in the 
approach, based on the use of braid statistics. They are characterised by the quantum di- 
mension, which is obtained as the solution of the mentioned discrete equations. We will give 
also some arguments in favor of stability of arising mappings of string nets, built of golden 
chain^. 



II. DIMER CONFIGURATIONS OF VORTEX PAIRS 



In this section we show how to construct dimer distributions starting from Hopf finks. 
Indeed, knots and links of field configurations appear naturally in the long-wavelength de- 
scription of {3+0)D systems with spin as follows. Let us consider a gauged Ginzburg-Landau 
model for the charged two-components order parameter $ = (0i,02)^ given by the free en- 
ergy functional 

F= [d'x\y— 

with a generic form of the potential V{(j)i, (j)2^ai.nd interacting with an internal vector gauge 
potential A. 

The components of the order parameter can be interpreted as the hopping and pairing 
amplitudes, respectively, in a t — J modeP^. The relations between these amplitudes can be 



2e 
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reshaped in terms of the function p and a real unit 3- vector n defined by 

= = (2) 

where a", a = 1, 2, 3 are the Pauli matrices. A residual phase-like degree of freedom is 
encoded into a momentum contribution to the effective gauge field 

2 

c = a-A, a = - ^ - c.c] . (3) 

Thus, one can map the model Q into the extended version of the 0(3) nonlinear a 
modeP^ 



(TX 



F = Fn + Fp + F, + = 
{dkuf + {dkpf + ^P^c^ + {Fik - Hikf + V{p, rii, ns) 



(4) 



Kaf3 = I d^xEikiaidkaf = j ^ | , where a] = at, = q. (6) 



where, in dimensionless variables, the field strength Fik = diCk — dkCi and the Mermin-Ho 
vorticity Hik = dittk — dkai = n ■ [din x dkn] have been introduced. 

Taking into account only homogeneous density states (i.e. p = const) free energy is 
bounded from below by the inequalitj^ 

+ + Fint ^ S2n' \Qf' (1 - \L\/\Q\Y , (5) 

given in terms of the entries of the symmetric matrixpHl 

Q L 
L Q' 

Taking the boundary condition n(oo) (0,0,1), one compactifies S"^ and Q is 

the degree of the mapping -n. : ^ S^. That is, Q measures the finking and knotting of 
the filaments, which are pre-images in of an arbitrary chosen value n{x,y, z) = hq E S"^. 
General arguments of the homotopic group theor}^ ensure that 7r3(S'^) = Z , i.e. Q e Z, 
labeUing distinct sectors among all configurations of the field n. For two finked loops one 
has (5 = 1, for the trefoil knot Q = Q and etc. On the other hand, L is the mutual linking 
numbei'^ of the fields c and a. The integer Q is relatecP^ to the coefficient k in the Chern- 
Simons action. For instance, the value Q = 1 for the Hopf link of two loops is equivalent to 
k = 2. 
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But, differently from n, the length of the vector c is not fixed and it does not belong 
to a compact manifold, thus homotopy group techniques say that L ( and Q') in ^ is an 
arbitrary real numbers. Furthermore, for \L\ < \Q\, from Eq. ([5]), a decreasing of the lower 



bound of F is obtained with respect to the case c = CP'. This effect is due to the interaction 
term —2FikHik in Eq. Q. Its contribution could be so important to completely compensate 
the energy contributions Fc and F^. Such a situation may occur because all contributions in 
Eq. ^ are of the same order. Then, one can prov^ the existence of inhomogeneous states 
with total energy smaller than in the Skyrme-Faddeev reduction (i.e. c = 0, p = const), if 
the knot is contained in a region of typical size > i? ~ 1/p and the amplitude of the 
momentum field |c| ~ cq is bounded by a/Ro < cq < I/-R0, where Rq can be figured out 
as the "thickness" of the filamentary structure of c, or n, and a = {Ro/R) < 1 is the 
packing parameter. Thus, by ^ one can evaluate the decrease of energy, by a negative 
(condensation energy) contribution AF ~ (647r^/(5^''^) a. 

In the optimum case of great values c the self-dual state with Fn = F^. may be considered 
as the ground state with Fmm = 0. It is characterized by the dense packing of filaments in 
knot^ and by the condition a — c = A = when the matrix 

i^«/3 = Q 1^ I j • (7) 

Indeed, let us assume that the amplitude of density p^ ~ R~^ is large enough and a 
knot size R is so small that the packing degree of filaments in the knolP^ a = ^/R < 1. 
The correlation lenght ^ having the order of the lattice constant is determined the filament 
thickness. The decrease of our tuning parameter a with the decrease of p^ is accompanied 
by the transition to a state, in which among planar field projections the field configurations 
in the form of closed one-dimensional distributions are most preferable. In accordance with 
Ginzburg's proposal, we will call this phase state with the spontaneous diamagnetism a 
toroid stat d^'^ ' ^^l The phase state with zero value of the total magnetic moment 
in which current order is characterized by a polar vector T, that changes sign under time 
reversal, can be represented by the ordering of toroid moments. The toroid moment T = 
I Xli ["^i ^ can be determined in the following way'^^'^ 

a = curl curl T . (8) 
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Using this equation and the identity!^ 

(curl a) . = ^ eiki n [dkn x din] (9) 
one can express the toroid moment via the degree of chirality n ■ [9fcn x diTL\ in the form: 



T,{v) = I dh' 

iOTT' 



r - r')k f ^3^„ n ■ [dkn x djii] 



(10) 



We see that toroid moment T is given by the Biot-Savart law as well as by the factor 
which is determined by the Coulomb Green function and by the field strenth Hik of the Hopf 
invariant density. The vector T characterizes the distribution of the poloidal component of 
the current on the toru^. The magnetic flux of this current is confined to the interior of 
the torus. If the degree n ■ [dkH x diu] = Ski 5(r) of the n-field noncoUinearity is locahzed 
at a point the toroid moment has the form Tj(r) = {^Tx)^^ekiXk/'r. The toroid moment is 
perpendicular to the plane, where the central loop of the torus is situated, and belongs to 
the perpendicular plane that intersects this loop. 

The significant constraint of the scale of toroid configuration appears from the analysis of 
stability of such field distributions. As it has been shown in Ref.^Sl, that these configurations 
are stable with respect to s-wave perturbations if the characteristic scale of linked distribu- 
tion is smaU enough, i.e. it is of the order of several lattice constants. Therefore, to describe 
topological states we have to use a lattice theory. Besides, in addition, the main contribu- 
tions of field configurations correspond usually to small values of topological numbers Q, 
L. Keeping all this in mind we can conclude that the Hopf links of loop pairs (see Fig. 1), 
that after being projected on the plane, provide dimer distributions of vortex pairs which 
are the the most significant configurations. To deal with states without the time inversion 
symmetry breaking, we should take into consideration totally neutral dimer configurations 
with antitoroidal ordering of the toroid moments. 
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Figure 1: (Color online). The Hopf pair with Q = 1. Their projection into the plane presents a 
system of two dimers, characterized by the respective orientation of toroid moments. 



III. DISCRETE EQUATIONS OF EXCLUSION STATISTICS 

Let us consider a system, which contains a set, {Na} of particles with types a. The 
collective index a = {a,i) contains the index a for denoting internal degrees of freedom 
and index i enumerates rapidities of particles. If we fix the variables of all particles, except 
the a-th one, the iV-particle wave function can be expressed via the one-particle function of 
the a-th particle. Let Da be the dimension of such a basis. Then the rate of changing the 
number of vacant states due to adding Ni, particles determined the matrix gab of statistical 
interaction in the following way 

Assuming that the matrix gab does not depend on the set of numbers {Na}, we have the 



solution of the Eq. (11): 



Da = -J29abNb + Dl (12) 
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The Eq. (12) contains the number of particles Nt^ added to the system, and the number of 
vacant states D° of the a-th type in the initial state without particles. The number of holes 
Da determines the statistical weight as follows 

{Na + Da-l + Eb9ab5abV. 



w 



n 



(13) 



{Nay.{Da-l + Eb9abSaby. ' 

In the cases gab = and gab = ^ab the Eq. ( p^ yields well-known statistical weights of Bose 
and Fermi particles. 

The statistical weight W allows to find the entropy S = \n.W and thermodynamical 
functions. The free energy in the equilibrium state 

F = -TY,DlHl + w-a^) (14) 

a 

is determined by the function Wa-, which can be found from the equation 

(1 + «^a) n (1 + ^^"') = e(^°-'^'')/^ (15) 

b 

The variable Wa can be expressed via so-called pseudo-energies e = T\n{Da/Na) by 



means of the parametrization Wa = e""/^. In Eq. (15), T, /!„ and e° are the temperature, 
the chemical potential and the bare energy of quasiparticles of the type a. 



We consider the solution of the equation (15) in the limit T e° — In this case the 



Eq. (15) can be written down in the form 

n(i+ 



') 



(16) 



which is typical for the thermodynamic Bethe ansatz. Here Nab = gab — Sab- 

Below we will be interested in the case of ideal statistic^S^^nZl, when phases of the scat- 
tering matrix, being functions of rapidities, have the structure of step functions. In this 



case the integral equation (16) transforms into the algebraic transcendental equation. The 
matrix Nab can be expressecpE^EHIil ^ia the incidence matrix Gab = 5a+i,6 + Sa,b+i of Lie 
algebra with the help of the identity = G{2 — G)^^. The matrix 2 — G is the Cartan 



matrix of the graph Ak+i/Z-^. Using this identity in Eq. (16) and replacing Wa = — 1 it 



is easy to see that the Eq. ( |16[ ) has the form 

1 + d2, a = 1, 



[k/2] 

di=i+ n 

j=l,b=2j 



1 + da-ida+1, a = 2, [k/2] - 1, 
1 + d\k/2\-id[k/2\, a = [k/2]. 



(17) 



Here index a is connected with the value of the spin j by the relation a = 2j, and the upper 
hmit of the product fixes the Jones-Wenzl projectoi'^^'^. We will show in Appendix that 
the Eq. (17) is presented in fact the special hmit of the Hirota equatiorP^. 



The distribution function 

_ 1 _ 1 _ 1 
'''^ " ^ " ~ e-/^ + 1 ^^^^ 

in our case coincides with the probability p{ad — > 0)''^ of annihilation of a particle-antiparticle 

pair in the system of two linked loops of world lines which describe the process of annihilation 

of two pairs. 



We can find the solution of the Eq. (17) taking into account the appropriate boundary 
conditions by comparing it with the identity 

[a]2-l = [a + l]Ja-l],. (19) 

Here [a\q = {q°- — q~°')/{q — q^^)-, q = e™/'''^^'^^ is the deformation parameter of the SU{2)k 



Chern-Simons theory. Identifying da = [a + l]g, we can see that solutions of the Eq. (17) 
are quantum dimension^^SE^E^ESlQ] 

_ sin[7r(a + l)/(fc + 2)] 

sin[7r/(fc + 2)] ' ^^^^ 

which are expressed via the Chebyshev polynomials of the second kind, Um = sin[(m + 

1)0]/ sinO with specification 9 = -K/{k + 2) for Ak+i algebra. In the limit A; 3> 1, da=2j equals 

2j ■ + 1. The meaning of the quantum dimension is as foUows. The quantum dimension da 

determines the rate d^ at which the dimension of the topological Hilbert space grows after 

particles are added. 

We pay our attention to the fact that the Eq. (17) is a fermionic representatioipJ^^^ 



of the recursion relation for the Chebyshev polynomials of the second kind. The bosonic 
representation of recursion relations has the form Um+i{x) + Um^i{x) — 2xUm{x) = 0. From 
this point of view, we can caU the Eq. ( p6| ) (see Appendix) the anyon representation of 
recursion relations for the Chebyshev polynomials of the second kind. The roots of the 
Chebyshev polynomials, being the eigen values of the matrix G, are equal to 

x^,. = g'"-^^ + g-^'-^^) = 2 cos (^^1^) • (21) 

The greatest eigen value xo^k of the incidence matrix G is given by the Baraha numbers 

d = 2cos['K/{k + 2)]. (22) 
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In particular for the special value = 3, we have the golden ratio d = {1 + v^)/2 which is 
the solution of the algebraic equation cP = d+ 1. For /c = 2, the Baraha number d and the 
quantum dimension d2j=i = V2 coincide. 

To clarify the meaning of the d's one should emphasize that (i) the numbers d determine 
eigen values of Wilson operator for the contractible unknotted loofP^. (ii) For special 
values of the parameter d = q + q~^, the generators Bi = I — qci satisfy the relations of 
a braid group under the condition, that the generators Cj satisfy the relation ef = dci of 
the Temperley-Lieb algebra, (iii) The values of the parameter d are nontrivial restriction, 
which leads to the finite-dimensional Hilbert spaces, (iv) The wavefunction \E', defined on 
the one-dimensional manifold, which is a joining up of the arbitrary tangle a and the Wilson 
loop Q, i.e. '^{aUQ), equals d'^{a'^. Thus the parameter d has the meaning of the weight 
of the contractible unknotted Wilson loop, and d"^ acquires the meaning of fugacitj^. (v) 
Besides, it turns out, that for the mentioned values of d, the theory is unitary. 

Summarizing one can say that the points of intersection of braid statistics and statistics 
with the generalized exclusion principle are the set of the Baraha points d = 2cos[7r/(A; + 2)], 
where the processes of braiding and fusion of string manifolds are self-consistently united. 



IV. DISCUSSION 

The problem of the coexistenc^ of locality and braiding can be solved by constructing 
Hamiltonians H = ^iHi of the Rokhsar-Kivelson (RK) typd^. Each term = QfQi 

( 1 

in the sum with Q = \ acts at the RK point as a projector builted by dimer 

\-i 1 ; 

configurations. If we locate the dimers on the opposite hnks of plaquetts we will encounter 
contradiction due to spatial separation of the braiding phenomena. The best way to solve 
this problem is the distribution of dimer states between odd and even sites of the latticd^. 
Another way corresponds to the use of the Fischer lattic^. The checkerboard distribution 
of linked dimer degrees of freedom with defects in the order of effective fluxes (Fig. 2) may 
be one of the possible ways to solve the problem. 

The reahstic candidate for the Hamiltonian with such a type of the ground state is the 
Hamiltonian Hi which contains the Temperley-Lieb generators in the form of the projectors 

= ^e, (23) 
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Figure 2: (Color online). A planar representation of the net. Shaded regions show distributions of 
the fluxes. 

to the singlet states. Obviously iff = Hi due to the Temperley-Lieb commutation relation 
= dci. The d's here are the Baraha numbers from the second section. Because of the rank- 
level symmetry, i.e. SU{N)k = SU{k)N, and the argument based on small values of integers, 
the SU (2)2 X SU (2)2 theory is a good candidate. The matrix of the 6j-symbol^in this case 



1 1 \ / 1 , 

and the braid operator is . Another important model 

i 



is equal to 

with computational universal rules is basecP^lon a golden chain builted by SU{2)^ x S'f/(2)3 
Fibonacci anyons. To construct 2D nets consisting of Fibonacci golden chains we should 
employ the Kirby calculus. This calculus is based on the application of Hopf hnks and is 
widely used in the theory of 3-manifolds. 

In summary, we found the quantum dimensions as exact solutions of discrete equations 
encoding brading and fusion processes. By means of projection of knotted field configurations 
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we proposed the quantum dimer models which incorporated braiding properties of one- 
dimensional topology. 
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V. APPENDIX 



Let us show that the Eq. (17) is a particular case of the Hirota equation^Sl 



T,%u + l)T^{u - 1) - T^,,{u)Tl,{u) = Tr\u)Tr'{u) 



(24) 



Here a is 
ra 



.v.iv. ^ the index of the agebra Ak+i, t is me uiscieie nme aim u is n. 
ipidities. The functions T^(u) are the eigen values of the transfer-matri?^ 
functions Y^^{u) = T^^-^{u)T^_-^{u) / [T^^^ {u)T^'^ {u)) the following equal: 

{l + Y,\,{u)){l + Y-_,{u)) 
iYr\u))-^){l + {Tr\u))-^ 



the discrete time and u is the discrete values of 

■ ■ For the gauged 



;ion 



^^'^ '^-(i + (r-V))-^)(i 

is vahd. In the Ai-algebra case the function Y^^{u) = Yt{u) satisfies the equation 



Ytiu + l)Ytiu - 1) = (1 + Yt+,{u))il + Yt^,{u)). 



(25) 



(26) 



Putting Yt = bf — 1 in this equation, in the hmit m ^ 1 we get the equation bf = 1 + bt+ibt^i 



which coincides with (17). We see that the function Yt equals Wt. 
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